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Abstract 

For a class of stationary regularly varying and weakly dependent multivariate time series 
(X n ), we prove the so-called complete convergence result for the space-time point processes 
of the form N n = M" =1 fi i / n .x,/a n )- As an application of our main theorem, we give a simple 
proof of the invariance principle for the corresponding partial maximum process. 
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1 Introduction 

Point processes theory is widely recognized as a useful and elegant tool for the extremal analysis 
of stochastic processes. This approach is splendidly illustrated by Resnick [IS] or Leadbetter and 
Rootzen [12]. It is well known for instance that for an iid sequence of random variables (V n ), 
regular variation of the marginal distribution is equivalent to the so called complete convergence 
result, that is to the convergence of point processes of the form 

n 

Nn — 52 ^(i/wAi/on)) ( 1 - 1 ) 

2=1 

towards a suitable Poisson point process. Such a statement then yields nearly all relevant asymp¬ 
totic distributional properties about the sequence (X n ), cf. [12]. There exists a rich literature 
on extensions of this result to dependent stationary sequences in both univariate and multivariate 
cases, for an illustration consider for instance Davis and Resnick [8], Davis and Hsing [6], Davis 
and Mikosch [7], Hsing and Leadbetter ra or Basrak et al. [3]. One of the earliest and key results 
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in this area was obtained by Mori [15] who showed that all possible limits for the point processes 
N n have the form of a Poisson cluster process provided that the random variables X' n s are strongly 
mixing with a strictly positive extremal index. In such a case the limit can be written as 


N - EE<Ww«>’ 

i j 


where JP ^(t ?:: p, ; ) is a suitable Poisson process and (^P §Q i:j )i represents an independent iid sequence 
of point processes on M\{0} with the property that max., \Qij\ = 1 for all i. However, the relation¬ 
ship between the sequence (X n ) and the distribution of the limiting process N or its components 
is up to now only partly understood. Due to the dependence, for a general sequence (X n ), there is 
no complete convergence result which determines the shape of the limiting process N in the form 
suggested by Mori. It is our main goal here to give such a result for a rather wide class of weakly 
dependent regularly varying processes, even in the multivariate setting. 

Recall that a d-dimensional random vector X is regularly varying with index a > 0 if there exists 
a random vector © G § d_1 such that 


1 

P(||X|| > x) 


P(||X|| > ux, x/\\x\\ g •) =>■ m"“P(0 g •), 


( 1 . 2 ) 


for every u > 0 as x —» oo, where denotes the weak convergence of measures. Note that the 
definition does not depend on the choice of the norm, i.e. if (II.2ft holds for some norm in M d , it 
holds for all norms, with different distributions of 0 clearly. A d-dimensional time series (X n ) n£Z is 
regularly varying if all of the finite-dimensional vectors ( X *.,..., Xi), k,l G Z are regularly varying, 
see Davis and Hsing [6] for instance. We will consider a strictly stationary regularly varying process 
(X n ) n£Z . This means in particular, that there exists a sequence (a n ), a n —> oo such that 


nP(||X 0 || > a n x ) —* x " for all x > 0. (1.3) 

According to Basrak and Segers [4], the regular variation of the stationary sequence ( X n ) is equiv¬ 
alent to the existence of the tail process {Y n ) n&z which satisfies P(||iKo|| > y) = y~ a for y > 1 and, 
as x — * oo, 

((oT 1 X n ) nez | \X 0 \ >x)^ ( Y n ) nez , (1.4) 

where —^0 denotes convergence of finite-dimensional distributions. Moreover, the so-called spectral 
tail process (© n ) ne z defined as a sequence G n = W n /||’K 0 || , n G Z, turns out to be independent of 
|| W 0 || and satisfies 

((HXoir 1 X n ) nez | ||Xo|| > x) ^ (0„) ne z. (1.5) 

as x — * oo. In the sequel, we assume that there exists a sequence (r n ), where r n — > oo and 
n/r n — > oo, such that (X n ) and (r n ) satisfy the following two conditions. Denote first O = 

I d \{0} = [—oo, oo]^ \ {0}. 


Main assumptions 
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{A!)\ For every / G C'^([0,1] x O), denoting k n = \n/r n \, as n —> oo, 


E 


exp 


2=1 


-n E 


k =1 


eX Pi 

i= 1 ' 


—>■ 0. 


( 1 . 6 ) 


(MC): For every u > 0, 


lim limsnpPI max ||Xb|| > a n u 

TTL fOO n —yoo \ 772^|2| n 


|Xn|| > Chr,u = 0 


(1.7) 


For thorough discussion of conditions {A!) and (AC) we refer to [3] or Bartkiewicz et al [I]. We 
observe here only that m-dependent multivariate sequences clearly satisfy both assumptions. More 
general strongly mixing sequences always satisfy condition A 1 , which roughly speaking, allows one 
to break the sequence into increasing and asymptotically independent blocks. The condition AC on 
the other hand restricts the clustering of extremes in the sequence (||.XL||) over time. In addition to 
m-dependent sequences, it is satisfied for many other frequently used heavy tailed models, including 
for instance stochastic volatility or ARCH and GARCH processes, see [I] and [3j. 

The main object of interest in the sequel is the multivariate version of the point process in (II.Tj) 

n 

N n — 'y ] $(i/n,Xi/a n ) 

2=1 

on the state space [0,1] x O. Throughout, the space of point measures on a given space, S say, is 
denoted by M p (S) and endowed by the vague topology, see Resnick [16]. 

Under very similar assumptions as above, Davis and Hsing in [6] analyzed the projection of the 
process N n on the spatial coordinate in the univariate case. They showed that the point processes 
K = E"=i Sxi/a„ j n G Z, on the space [— oo, 0) U (0, oo], converge in distribution to a point process 
which has a Poisson cluster structure which is only implicitly described. They further apply this 
result to determine the limiting distribution of the partial sums in the sequence (X n ). These results 
were extended to the multivariate setting by Davis and Mikosch [7j. Clearly, unlike N n , the point 
processes N* contains no information about time clustering of extremes in the sequence (X n ). 

Point processes N n were already studied by Basrak et al. in [3] in the univariate case, but the proof 
therein carries over to the multivariate case as observed in Basrak and Krizmanic [2] with some 
straightforward adjustments, to show 


AL 


[0,l]x(R\[—u,u]) d 




[0,1]x(R\[— 


( 1 . 8 ) 


as n —> oo, for any threshold u > 0, with the limit N bd unfortunately depending on that threshold. 
This asymptotic results can be used to deduce functional limit theorems for partial sums or partial 
maxima in the time series (AL), see Theorem 3.4 in [3] or Proposition 14.11 below. Our goal here 
is to avoid restriction to various domains in (1 1.8 j) . to find the correspondence of this result with 
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the earlier results of Davis and Hsing [6] and Davis and Mikosch [7], Our main theorem below 
essentially unifies theorems on the limiting behavior of point processes given in [HI S]. Moreover, 
it reconciles their apparently very different statements in the framework suggested by Mori’s result 

m- 

In the following section we show an interesting preliminary result about the structure of extreme 
clusters in the process (X n ). Then, in Section [3] we prove a general theorem about point process 
convergence for regularly varying time series. In Section [4] we apply this theorem to show the 
invariance principle for the so-called maximal process in the space D[0,1]. A corresponding theorem 
for iid sequences is well known and can be found in Resnick |l6j . For nonnegative stationary 
regularly varying sequences, it was recently proved by Krizmanic m ■ Our version of this result 
includes other regularly varying sequences, and since it relics on our main theorem, the proof is 
straightforward and relatively simple. We also exhibit a simple technique to avoid problems at the 
left tail which are typically a source of frustration in this sort of limiting theorems. 


2 Preliminaries 

It was shown in [4] that under the main assumptions above, there exists 


9 = P sup ||10|| < 1 = lim lim P (M r < x ||X 0 || > x) > 0, 

' / 1 —>-oo x —^oo 


( 2 . 1 ) 


where M r = max{||AA|| : k — i,... , r}. The number 9 > 0 represents the extremal index of the 
random nonnegative sequence (||X t ||). 

Following |4], we define an auxiliary sequence (Zj)j G % as a sequence of random variables distributed 
as (Y j)j e z conditionally on the event {sup i< _ 1 ||1A|| < 1}. More precisely, 


c (E s z) = c ('E s ' 

\j£ Z / VieZ 


sup ||Wi|| < 1 I , 

i <-1 


( 2 , 2 ) 


where jC(N) denotes the distribution of the point process N. It was shown in Theorem 4.3. in {4] 
that 


^1 

' i=l 


Afr n a n u ) ^ ^ &Zj ■ 

j&Z 


(2.3) 


Expressing the weak convergence of point processes using Laplace functionals, (12. 3 p is equivalent 
to 

(2.4) 


L-HZlf^nU) 'Xi) I M > a A _ > E 

e f(Yj) 1 SU p H^ll < i 

V 1 / 

j<-i 


for all u G (0, oo) and / G C 


I< 
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Observe that by AC, ||W„|| —* 0 with probability 1 as \n\ —> oo, see Proposition 4.2 in [4|. Therefore 
the same holds for Z n J s in (12.21) . In particular, the random variable 


L z = sup \\Zj\\ 

jez 

is a.s. finite, and clearly not smaller than 1. To determine the distribution of L z , observe that 
from Proposition 4.2 in [4] it follows that 

k n ¥(M rn > a n u ) —» 9u~ a . (2.5) 


Now for v > 1, we have by (I2.3j) 


¥(Lz >v) = l -P J oo)) = 0 


1 - lim P y'<5(a tlU )-H|A: i ||((v, °°)) = 0 I M r n > a n u 


, 2=1 


i , r . /c„P(M, > a n uv) 

= lim P (M rn > a n uv \ M rn > a n u ) = lim 


n^oo knP(M rn > a n u) 


Therefore, 


P(L Z > v) — v ". 


For ( Zj ) as in (12.2p . we dehne a new sequence (Qj)j^z by 


( 2 . 6 ) 


Qj — Zj/L z , j G Z . 

We will show the independence between the point process and the random variable L z , 

which might not be entirely surprising in the view of the independence between ||Wo|| and the 
spectral tail process in (11.51) . However, this result seems to be new. 


Proposition 2.1 Assume that a regularly varying stationary sequence ( X n ) with corresponding 
sequence (r n ) satisfies conditions A' and AC. Then 


, i =1 


M r . 


MpfXii 


a„u 


Mr n > ci n u 



Moreover, L z and on ^e right hand side are independent. 


(2.7) 


Proof. For m = JT S Xi G M p (0) denote by x m the largest norm of any point in m, i.e. 
x rn = supj ||*j|| and dehne the mapping f by 

4> : m hg {rn, x m ). 
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As observed in [6] such a mapping is continuous. Hence, applying 0 to (12.30 . we obtain 




M r 




. i=1 


a^u 


Mr > a„u 



Consider for v G M p and b G (0, oo), the mapping 

if : {y, b) u b G M p , 


( 2 . 8 ) 


where z/j,(-) = v(b 1 -). Mapping if is again continuous by Proposition 3.18 in Resnick [T7] for 
instance. Hence, (12.8p implies (12.7j) by the continuous mapping theorem. 

To show the independence between Lz and )T7 Sq , it suffices to show 


E 


exp \-^2fi{Qj) ) 1 


{L z >v} 


= E 


exp ( _ 5Z A(Qj 


> u) , 


(2.9) 


for an arbitrary function f\ G (7^(0) and v > 1. By (12.711 . the left-hand side of ({279]) is the limit of 


E 


eXP ( “ f'( X i/ M rn) ) 1 


{(ua n )- 1 M rn >v} 


2=1 


A/j- > a n 


which further equals 


E 


exp ( fi(Xi/M rn ) 


2=1 


M r „ > a n v 


P(M r „ > q n u) 
P(M rn > a n ) ' 


By (12.71) and the continuous mapping theorem, the first term above tends to E exp h(Qj) 

as n —> oo. By (12.5ft . the second term tends to P (Lz > v) = v~ a . which implies (12.9ft . 


3 Main theorem 

Theorem 3.1 Let ( X n ) ne z be a stationary series of jointly regularly varying random vectors with 
index a, satisfying conditions A! and AC for a certain sequence ( r n ). Then 

n 

N n = ^ S(i/n, Xi /a n ) ^ N = ( 3 ‘ 1 ) 

*= 1 * 3 


where 

i) Yhi d(Ti,Pi) a Poisson process on [0,1] x (0, oo) with intensity measure Leb x v where u(dy) = 
Qay~ OL ^ 1 dy for y > 0. 


6 





























ii) (Ejtvij)i an i-i-d- sequence of point processes in 0 independent of )T /l b{T t ,p,) and with 
common distribution equal to the distribution of Yhj ■ 

Proof. To show (14.2p . it is sufficient to show that 

Ee -^(/) ^ Ee- N{f) , 

for an arbitrary / G C^QO, 1] x O). Observe that for any such function / there is a constant u > 0, 
such that the support of / is a subset of [0,1] x O u . By (jl.8p . i.e. by the multivariate expansion of 
Theorem 2.3 in [3], we know that 


Ee-^i —> 


exp 


- / (1 - Ee " f ^ uZ ^) du~ a dt 


(3.2) 


Consequently, it suffices to show that the right hand side above corresponds to Ee for any 
such function / and corresponding u > 0. 

As in Theorem 2.3 in [3], one can write 
Ee -mf) = Eexp 

* 3 

E ( P exp{- V f(Ti, PtViMin, Pt))> 

\ i j 

E ( exp{- V f(T it PivPMn, Pt)), 


= E 


= E 


where the last equation follows from Lemma 3.10 in [16], since Yhi^{Ti,Pi) and on the 

right-hand side of (14.2ft are independent. Define h(t,v) = Eexp{— Y2j f(t,vQj)}. We have 


Ee _JV(/) = Eexp ^ In h(T u Pf 


The right-hand side is the Laplace functional of )Tb S^pa evaluated at — In h. Since JT 5( Ti,Pi ) ~PRM(Le6x 
u) on [0,1] x (0, oo), the Laplace functional of this process has a very special form, see Resnick [16] . 


Ee-*W = 


exp 


= exp 


/ (1 — h(t, v)) (Leb x v)(dt, dv) 

qo,i]x(o,oo) 

a oo 

(1 - Ee~^ f(t ’ vQ ^) dav~ a ~ 1 dvdt 

. 


(3.3) 
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Consider now the right hand side in (13.21) . Using Z :j = L Z Q p the distribution of L z calculated in 
(12.61) and independence shown in Proposition 12.11 we have 


fi 


exp 

= exp 

= exp 

= exp 


1 POO 



- I (1 - Ee“ f(t ’ uZ ^) 9u~ a dt 

(1 - E al~ ot ~ 1 dWu~ a dt 
(1 - Ee-^ f(t ’ vQ >y) av-^dvddt 
(1 - Ee _E 3 f d, v Qj)') 6av~ a ~ 1 dvdt 


'0 J 1 

poo 



' 0 J u 

poo 



0 JO 


= Ee _JV(/) , 


where the second equality follows using the change of variable v = ul, and the last equality follows 
from the fact that sup^ ||Q. ; || = 1 and f(t,x ) = 0 for x < u. ■ 


As we discussed in the introduction, one consequence of the theorem above is the functional limit 
theorem for partial sums of the univariate sequence (Ah), see Theorem 3.4 in [3]. In the multivariate 
case, under the conditions of Theorem 13. II it was shown by Davis and Mikosch [T] that partial sums 
S n = X 1 + --- + x n , n > 1, satisfy 


— 

dn 


OL 1 


(3.4) 


for a G (0,1) and an a-stable random vector £ Q , actually the same holds for a G [1, 2) if one assumes 
for instance that X^s have a symmetric distribution and that for any S > 0 the following standard 
technical assumption holds lim^olimsup^gQPdla” 1 X i I{||x,||<£a n }|| > S) = 0, see [7] again. 


As observed by Mikosch and Wintenberger [13] , under condition that 



Q 


< oo, 


which indeed always holds in the case a < 1 or if (X n ) is m-dependent, the spectral measure T Q 
of the stable random vector in (13. 4 j) can be characterized as follows 



where £ G S d x , (t, s) denotes the inner product on on S d 1 and || • || denotes the Euclidean norm. 
For a definition of a the spectral measure of a stable random vector, see [ 19l Section 2.3]. This 
can be used in certain cases to give an alternative representation of the so-called cluster index 






















b(t) G S> d 1 studied by Mikosch and Wintenberger in [131 H3] • Indeed, for a regularly varying (X n ) 
which is Markov chain under conditions of [IT, Theorem 4.1] 


b{t) 


1 — a 

T(2 — a) cos(7ra/2) 



{t, s)°lT a (ds). 


As observed in mm, the cluster index b(t) plays the key role in the description of large deviations 
for partial sums of regularly varying sequences. 

A more direct consequence of Theorem 13.11 is the functional limit theorem for the partial maxima 
of a univariate sequence considered in the following section. 


4 Invariance principle for maximal process 

Denote by M' n = max{Ah,..., X n }, n > 1. Following Resnick [16] or Embrechts et ah [9] we 
consider the following continuous time cadlag process 

V m = / M '[nt\/ a n t>l/n, 

n[) \ x 1 /a n t<l/n, 

indexed over the segment [0,1]. It is actually customary to exclude t — 0 to avoid technical 
problems, see [TB]. In our approach, such issues are easily avoided, and therefore we include t — 0. 
From the practical perspective, it seems sufficient to consider t < 1. Extending the main result 
below to t G [0, oo) remains a technical, but relatively straightforward task, see Chapter 3 in 
Billingsley [5]. 

Recall that the extremal process generated by an extreme value distribution function G (G-extremal 
process, for short) is a continuous time stochastic process with finite dimensional distributions 
G si ,...,s k satisfying 

G n ,...,s h (x i, ...,x k ) = G S1 (/\i =1 Xi)G S2 - S1 (A*L 2 £j) • ■ (x k ) , 

for all choices of k > 1, 0 < si < • • • < s k , Xi G M, i = 1,..., k, see Resnick [TB], The processes Y n 
introduced above are clearly random elements in the space of real valued cadlag functions D[ 0,1]. 
We will show that they converge weakly to a particular extremal process with respect to Skorohod’s 
Mi topology on D[ 0,1]. We refer to Whitt [20] for the dehnition and discussion of various topologies 
in that space, see also [3|. 

In the proposition below we make a small technical assumption about the right tail of the marginal 
distribution of X t ’s, i.e. we suppose that it satisfies liminf n _ ) . oo nP(X 0 > a n ) > 0. This implies that 

nP(A 7 ' 0 > a n x) —> px~ a for all x > 0, (4.1) 


and some constant p G (0,1]. 
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Proposition 4.1 Let (X n ) ne % be a stationary sequence of jointly regularly varying random vari¬ 
ables with index a, satisfying conditions A' and AC. Assume that liminf^oo nP(X 0 > a n ) > 0, then 

Y n =► (4.2) 

where £(£), t > 0 is a G-extremal process for the nonstandard Frechet distribution function G{x) = 
e~ KX “, x > 0, for some n > 0, and the convergence takes place in D[0,1] endowed with the 
Skorohod’s Mi topology. 

Proof. Consider the functional T + : M p ([0,1] x O) —> Z?[0,1] on the space of Radon point 
measures given by 

T + (m)(t ) = sup ji V 0 , (4.3) 

ti<t 

where m = JR denotes an arbitrary Radon point measure on the space [0,1] x E, where we 
set for convenience sup 0 = 0. 

Denote M' p = {m £ M p : m([0, s] x (0, oo]) > 0 for all s > 0} and M” = {m £ M v : m({0,1} x 
(0, oo]) = 0}. We will show that T + is continuous on the set M p DM p . Assume, m n A m £ M' p C\M p . 
Because, m is a Radon point measure, the set of times t for which m({t} x O) = 0 is dense in [0,1], 
For all such Cs 

T + m n (t) —> T + m(t ), 

as n —> oo by Proposition 3.13 in pR|- Observe that T + m is a nondecreasing function for any m, 
therefore an application of Corollary 12.5.1 in Whitt |20| yields the convergence of T + m n —> T + m 
in D[0,1] endowed with Mi topology. 

Observe now that the limiting point process N in Theorem 13.II lies in M p fl M p a.s. Hence for N n in 
the same theorem, an application of the continuous mapping theorem yields the weak convergence 

in 

T + N n => T + N , 

as n —> oo. 

By Slutsky argument, to show that Y n =>• T + A^, it is sufficient to prove that 

d Ml (Y n , T + N n ) —> 0 , 

where djvri denotes the Mi metric on the space £>[0,1], see Whitt [20] . This follows from the fact 
that metric d.Mi is weaker than the uniform metric d^, and the following obvious limit 

doo{Y n , T + N n ) = o} "A 0 , 

a n 

as n —> oo. 

Observe that, by definition, iV([0, t] x (0, oo]) = 0 implies T + N(t ) = 0 for a fixed t > 0. Therefore 

T + N(t ) = sup Pi sup rjij = sup PiJJi , (4.4) 

Ti<t 
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denoting U t = sup,,- ryj V 0. Note further that U-s form an iid sequence of random variables on the 
interval [0,1]. By (14. ip . P(f/j G (0,1]) > 0. Moreover, sequence (U t ) is independent of the PRM 
Sj £ Ti,Pi■ It is straightforward to check that JT £r i ,p i t/ i is PRM with mean measure Leb x z/ where 
n'(dy) = 9¥JJ a ay~ a ^ 1 dy for y > 0, cf. propositions 3.7 and 3.8 in Resnick [16]. It follows, by 
proposition 5.4.4 in Embrechts et al. [9] for instance, that T + N(t), t > 0 in (14.41) is a G-extremal 
process for 


G(u) = P (T + N( 1) < u) 


P 



exp (-em a u ~ a ), 


for any u > 0. 


By the proof, the constant n in Proposition 14.11 equals 

k = em a , 

where U = sup jTjij V 0. In the iid case, the extremal index 9 = 1, while JN 5 Vij = Sq where 
P (Q = 1) = 1 — P(Q = —1) = p. Therefore, k = p as known from Proposition 4.20 in Resnick [16] . 
In the case of nonnegative random variables p = 1 and sup ryj = 1 a.s. Therefore k = 9 in this 
case, cf. Krizmanic HU- 

Remark 1 Observe that more commonly used J\ topology is not applicable in our setting, because 
of the clustering of extremes. For an illustration of the problem, consider the point measures 
rn n = < 5 i/ 2 —i/n,i/2 + 81 / 2,1 -» m = hi/2,1/2 + £1/2,1 for n ->■ 00, and note that for T + in ( 14 . 3 ) 1 , T + m n 
does not converge to T + m in J\ topology. 

In the case 9 = 1, the limiting point process N is a Poisson random measure, and the convergence 
in the theorem above holds in the standard .J\ topology. The proof only has to be adapted to show 
that T + is an a.s. continuous functional with respect to the distribution of such a process. Such a 
result was already stated in Remark 2 of Mori m- 

Acknowledgements: This work has been supported in part by Croatian Science Foundation 
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